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Abstract 

The current observations are usually explained by an accelerating ex- 
pansion of the present universe. However, with the present quality of 
the supernovae la data, the allowed parameter space is wide enough 
to accommodate the decelerating models as well. This is shown by 
considering a particular example of the dark energy equation-of-state 
Wfj, = p<f>/ Pcf> = —1/3, which is equivalent to modifying the geometrical 
curvature index k of the standard cosmology by shifting it to (k — a) 
where a is a constant. The resulting decelerating model is consistent 
with the recent CMB observations made by WMAP, as well as, with 
the high redshift supernovae la data including SN 1997ff at z = 1.755. 
It is also consistent with the newly discovered supernovae SN 2002dc 
at z = 0.475 and SN 2002dd at z = 0.95 which have a general tendency 
to improve the fit. 

Key words: cosmology: theory, dark energy, CMB observations, SNe 
la observations. 

1 INTRODUCTION 

It is generally believed that the expansion of the present universe is acceler- 
ating, fuelled by some hypothetical source with negative pressure collectively 
known as 'dark energy '. This belief is mainly motivated by the high redshift 
supernovae (SNe) la observations, which cannot be explained by the decel- 
erating Einstein deSitter model which used to be the favoured model before 
these observations were made a few years ago. (It may, however, be noted 
that if one takes into account the absorption of light by the inter-galactic 
metallic dust which extinguishes radiation travelling over long distances, then 
the observed faintness of the extra-galactic SNe la can be explained success- 
fully in the framework of the Einstein deSitter model. This issue will be 
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discussed in a later section. In the rest of the discussion, we shall not include 
this effect.) 

Although the best-fitting standard model to the SNe la data predicts 
an accelerating expansion, however, it should be noted that the low density 
open models, which predict a decelerating expansion, also fit the SNe la 
observations reasonably well. Unfortunately, these models are ruled out by 
the recent measurements of the angular power fluctuations of CMB including 
the first year observations made by WMAP. 

In this paper we show, by considering a particular example, that both 
these observations (and many others too) can still be explained by a decel- 
erating model of the universe in the mainstream cosmology. 

The observed SNe la explosions look fainter than their luminosity in the 
Einstein-deSitter model. This observed faintness is explained by invoking a 
positive cosmological constant A following the fact that the luminosity dis- 
tance of an object can be increased by incorporating a 'matter' with negative 
pressure in Einstein's equations. However, a constant A, is plagued with the 
so called the cosmological constant problem: Why don't we see the large 
vacuum energy density p v = A/8nG = pp\ 10 76 GeV 4 , expected from par- 
ticle physics which is ^10 123 times larger than the value required by the SNe 
observations? Or, why are the radiation energy density and the energy den- 
sity in A (A/87rG) set to an accuracy of better than one part in 10 123 at the 
Planck time in order to ensure that the densities in matter and A become 
comparable at precisely the present epoch? (Padmanabhan 2002; Peebles 
& Ratra 2003; Sahni & Starobinsky 2000). A phenomenological solution to 
understand the smallness of A is supplied by a dynamically decaying A. A 
popular candidate is an evolving large-scale scalar field 0, commonly known 
as quintessence, which does not interact with matter (T^ ■ = 0) and can 
produce negative pressure for a potential energy-dominated field (Peebles 
& Ratra 2003). This is equivalent to generalizing the equation of state of 
vacuum to 

P<t> = (1) 

Moreover, it is always possible to add a term like on the right hand side 
of Einstein's field equations independent of the geometry of the universe. 
The equation-of-state parameter is, in general, some function of time and 
leads to a constant A for w$ — — 1. 

It would be worthwhile to mention here that there is another approach 
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to understand the smallness of A by invoking some phenomenological models 
of kinematical A which result either from some symmetry principle (Vish- 
wakarma 2001a, 2001b), or from dimensional analysis (Vishwakarma 2000, 
2002a; Carvalho, Lima & Waga 1992; Chen & Wu 1990) or just by assum- 
ing A as a function of the cosmic time t or the scale factor S(t) of the 
R-W metric (for a list of this type of models, see Overduin & Cooperstock 
1998). Although the quintessence fields also behave like a dynamical A (with 
Aeffective = SnGp^), they are in general fundamentally different from the 
kinematical A models. In the former case, quintessence and matter fields 
are assumed to be conserved separately, whereas in the latter case, the con- 
served quantity is [T^ atter + (A(t) /87rG}(? u ], which follows from the vanishing 
divergence of the Einstein tensor. This characteristic of the kinematical A 
models makes them consistent with Mach's principle (Vishwakarma 2002b). 
However, in the following, we shall keep ourselves limited to the former case 
with a constant w^. 

2 FIELD EQUATIONS AND DISTANCE MEASURES 

In order to study the constraints on the parameters coming from the different 
observations, we shall describe the model briefly in the following. For the 
R-W metric, the Einstein field equations, taken together with conservation 
of the energy, yield 



H\z) = H t 



an0(l + ^) 3 + ^r0(l+^) 4 + ^0(l+^) 3{1+ ^ ) -^0(l + ^ 



(2) 



q{z) 



HI 
H 2 (z) 



9 - m0 'i + zf + a (i + zY + (1+ o 3w ^ (i + ^) 3(1+ ^ 



2 v 2 

'(3) 

where fL £1X6 . clS usual, the different energy density components in units of 
the critical density p c = 3H 2 /8nG: 

n __ 8nG = 8nG = 8nG 

~ 3iJ 2 ^ m ' — 3i7 2 ^ r ' — 3iJ 2 ^' — S 2 H 2 

and the subscript denotes the value of the quantity at the present epoch. 
Equation (2) at z = implies that Q m0 + f2 r0 + fi^o = 1 + ^ko, by the use of 
which Equation (2) reduces to 
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H\z) = Hi ^{l + zf + ^{l + zY + ^{l + z)^ 1+w ^ 



+ (1 — ^mO ~~ ^rO — Q<t>o) (1 + 2 



(5) 



In a homogeneous and isotropic universe, the different distance measures 
of a light source of redshift z located at a radial coordinate distance r\ are 
given by the following. 
The luminosity distance d^: 



d L = (1 + z)S n, 
The angular diameter distance dp^: 



d\ = 



Son 
(1 + *)' 



(6) 



(7) 



where r\ is given by 



sin (i /o ^)) , when fc = 1 

i/o#y, when k = 
[ ™ h (i/o 2 ]|l). when k = -l. 



(8) 



The present value of the scale factor So, appearing in equations (6 - 8) which 
measures the curvature of spacetime, can be calculated from equation (4) as 



So 



Ho\ 



k 



\ (^mO + ^rO + ^00 — 1) 



(9) 



3 DEGENERACIES IN FRIEDMANN MODELS COMING FROM 
SNe la OBSERVATIONS 

Let us start with the SNe la observations which directly measure the ap- 
parent magnitude m and the redshift z of the supernovae. The theoretically 
predicted (apparent) magnitude m(z) of a light source at a given redshift z 



4 



is given in terms of its absolute luminosity M and the luminosity distance 
o?l through the relation 



m(z)=log w [D L (z)]+M, (10) 

where £>l = -£Wl is the dimensionless luminosity distance and M. = M — 
5 log 10 if + constant (the value of the constant depends on the chosen units 
in which d L and H are measured. For example, if d L is measured in Mpc 
and H in km s -1 Mpc -1 , then this constant comes out as m25). 

We consider the data on the redshift and magnitude of a sample of 54 SNe 
la as considered by Perlmutter et al (excluding 6 outliers from the full sample 
of 60 SNe) (Perlmutter et al 1999), together with SN 1997ff at z = 1.755 
(m eS = 26.02 ± 0.34), the highest redshift supernova observed so far (Riess 
et al 2001; Narciso et al 2002). In addition to this old sample of 55 SNe, we 
also consider the two newly discovered supernovae SN 2002dc at z — 0.475 
(m cff = 22.73 ± 0.23) and SN 2002dd at z = 0.95 (m cff = 24.68 ± 0.2), which 
have been discovered with the help of the recently installed Advanced Camera 
for Surveys on the Hubble Space Telescope (HST) (Blakeslee et al. 2003). In 
the following, we shall also keep an eye on how the addition of these SNe to 
the older sample affects the fit. 

The x 2 value is calculated from its usual definition 



N 



x 2 = E 



m eS _ m ( Zi ) 



(11) 



i=i 

where N is the number of data points (SNe) considered, m'f refers to the 
effective magnitude of the ith SN which has been corrected for the SN 
lightcurve width-luminosity relation, galactic extinction and K-correction. 
The dispersion 5m^ s is the uncertainty in m| ff . 

If one minimizes \ 2 in the standard cosmology (w^ = —1) by varying 
the free parameters fi m o, fi^o — ^ao an d M. (the contribution from the term 
containing fi r0 is negligible for this dataset and can be safely neglected), 
one finds that the best-fitting solution, from the older sample of 55 SNe, is 
tt m0 = 0.67, fi A o = 1-24 and M = 23.92 with x 2 = 56.93 at 52 degrees of 
freedom (dof) (= number of data points iV— number of fitted parameters), 
i.e., x 2 /dof = 56.93/52 = 1.09, a good fit indeed. Addition of the new 
SNe to the older sample gives a similar best-fitting model: f2 m0 = 0.64, 
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fi A0 = 1-20 and M = 23.92 with x 2 /dof = 57.78/54 = 1.07, showing a slight 
improvement in the fit. 

A 'rule of thumb ' for a moderately good fit is that x 2 should be roughly 
equal to the number of dof. A more quantitative measure for the goodness- of- 
fit is given by the x 2 -probability which is very often met with in the literature 
and its compliment is usually known as the significance level (should not be 
confused with the confidence regions). If the fitted model provides a typical 
value of x 2 as x at n dof, this probability is given by 

I roc 

Q ^ = VUVA , e ~ Uun/2 ~ ldu - (12) 
1 (n/2) Jx/2 

Roughly speaking, it measures the probability that the model does describe 
the data and any discrepancies are mere fluctuations which could have arisen 
by chance. To be more precise, Q(x,n) gives the probability that a model 
which does fit the data at n dof, would give a value of x 2 as large or larger 
than x. If Q is very small, the apparent discrepancies are unlikely to be 
chance fluctuations and the model is ruled out. 

The probability Q for the best-fitting standard model comes out as 29.7% 
(without including the new points) and 33.7% (with the new points), which 
represent good fits. These models represent accelerating expansion of the 
present universe according to equation (3), which suggests that go < if 

< —1/3 and fl^o > f2 m0 /(3 | | —1). Thus the best-fitting flat model 
(fl m + fl\ = 1) is also an accelerating one, which is obtained as Q m0 = 
1-Qao = 0.31 with x 2 /dof = 58.97/53 = 1.11 and Q = 26.6%, from the older 
sample of 55 SNe. Addition of the new points improves the fit marginally by 
giving n m0 = 0.32 with xV dof = 59.67/55 = 1.08 and Q = 31%, as the new 
best-fitting solution. 

Let us now examine the status of the decelerating models, particularly, the 
models with = 0. An important model in this category is the canonical 
Einstein-deSitter model (fi m o = 1, ^ao — 0), which has a poor fit: x 2 /dof 
= 93.01/54 = 1.72 with Q = 0.08%, from the older sample of 55 points. 
Even the addition of the new points does not improve the fit significantly, 
giving x 2 /dof = 94.61/56 = 1.69 with Q = 0.1% and the model is ruled out. 
However, the open models with low Q ni0 have reasonable fits. For example, 
the following decelerating models with A = (even without including the 
new points) provide: 

Ono = 0.2: xVdof = 66.74/54 = 1.24, Q = 11.4%; 
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n m0 = 0.3: xVdof = 68.78/54 = 1.27, Q = 8.5%; 
tt m0 = 0.4: xVdof = 71.31/54 = 1.32, Q = 5.7%, etc., 

which have reasonable fits and by no means are rejectable. Addition of the 
new points improves the fits marginally, giving: 
n m0 = 0.2: xVdof = 67.04/56 = 1.2, Q = 14.8%; 
il m0 = 0.3: x 2 /dof = 69.09/56 = 1.23, Q = 11.2%; 
tt m0 = 0.4: xVdof = 71.69/56 = 1.28, Q = 7.7%, etc. 

Open decelerating models with non-zero A (0 < Qao < ^mo/2) have even 
better fit. The best-fitting model with a vanishing A is obtained as Q m o = 
with x 2 = 64.6 at 54 dof (x 2 /dof = 1.2) and Q = 15.3%, from the older 
sample of 55 points. Addition of the new points to this sample improves the 
fit by giving x 2 /dof = 65.15/56 = 1.16, Q = 18.8%. 

Figure 1 shows the allowed regions by the full data in the f2 m0 — M. plane 
at different confidence levels for a vanishing A cosmology. Thus the low f2 m0 
open models with a vanishing A are fully consistent with the current SNe la 
observations. This result is also consistent with the findings of Gott et al 
(2001) from their median statistics analysis that the open models with low 
density and vanishing A are not inconsistent with the SNe la data. Addi- 
tionally, a low Q m0 is also consistent with the recent 2DF and Sloan surveys 
which give Q m0 = 0.23 ± 0.09 (Hawkins et al. 2002) and fi m0 « O.Uto.H (at 
2a) (Dodelson et al. 2001). However, these models do not seem consistent 
with the first-year observations of the temperature angular power spectrum 
of the CMB, measured accurately by NASA's explorer mission "Wilkinson 
Microwave Anisotropy Probe" (WMAP) (Bennett et al., 2003), or even with 
the earlier CMB observations (de Bernardis et al. 2000, 2002; Lee et al. 
2001; Halverson et al. 2002; Siever et al. 2002). These observations, when 
fitted to the standard cosmology, appear to indicate that f2 m0 + Qao ~ 1. 

4 THE PROPOSED MODEL 
4.1 Motivation 

We shall now consider a particular model of dark energy specified by the 
equation of state = —1/3, which, as we shall show in the following, 
explains both the observations - SNe la, as well as, CMB - very well. The 
equation of state = —1/3, for which varies as 1/S 2 (say, = a/S 2 H 2 ), 
is interesting in its own right. For this case, the Hubble and the deceleration 
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Figure 1: The allowed regions by the SNe la data (with 57 points) are shown in 
the O m o — M plane for = cosmology (A4 is usually termed as the 'Hubble 
constant-free absolute luminosity 1 ). The ellipses, in the order of increasing size, 
correspond to respectively 68%, 95% and 99% confidence levels. 
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(13) 



(14) 



which describe a decelerating expansion. Interestingly in this case, fi ma tter 
(= Q m + Q T ) alone is sufficient to describe H(z) completely, which gives the 
expansion dynamics of the universe. The additional knowledge of f^ decides 
the curvature through fl k0 = Q mQ + Q r0 + Q^q — 1. Expressions (13) and (14) 
are the same as the ones in the standard cosmology with A = 0, with only 
one exception: in the standard cosmology, f2 m0 + fi r o = 1 + &ko, whereas in 
the present model, Q m0 + Q r0 + Q^o = 1 + Qko- As mentioned earlier, f^ 
contributes to the curvature through So = Hq 1 J k / (ytmQ + fi r0 + ^</>o — 1) 
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and hence to the different distance measures for k — ±1 cases (for k — 0, 
distances do not depend on So). Although f2^ does not contribute to the 
expansion dynamics of the model (and neither to the distances for k = 
case), it helps Q m0 to assume even those values which are not allowed in 
the standard cosmology, and hence gives more leverage to fi m0 . It is in fact 
this property of the model which is instrumental in explaining the CMB 
observations even for low fl m o, as we shall see in the following. 

It may be noted that although the topological defects, like cosmic strings 
and textures, also have an equation of state p + 3p = (i.e., their density falls 
off as 1/S 2 ), however, the converse of this is not true and a dark energy with 
Wfj) = —1/3 need not necessarily be represented by cosmic strings. Moreover, 
this 'pseudo-source' term does not explicitly contribute to the expansion 
dynamics (and hence to the deceleration or to the expansion age) of the 
model but essentially to the curvature, as we have mentioned earlier. One 
can obtain the same cosmology by removing this term and simultaneously 
modifying the curvature index k by {k — a) (Vishwakarma & Singh 2002). 
It would, therefore, be more appropriate to consider it as a shift in the 
geometrical curvature of the standard cosmology and not as a source term. 
This term also appears in the context of brane cosmology by adding a surface 
term of brane curvature scalar in the action (Vishwakarma & Singh 2003; 
Singh, Vishwakarma & Dadhich 2002). 

4.2 CMB Observations 

The temperature fluctuations AT in the angular power spectra of the CMB 
correspond to oscillations in the AT - £ (Legendre multipole) space. The 
peaks of these oscillations can be explained in terms of the angle 9a sub- 
tended by the sound horizon at the last scattering epoch when CMB photons 
decoupled from baryons at Zd ec = 1100 (Hu & Dodelson 2002). The power 
spectrum of CMB can be characterized by the positions of the peaks and 
ratios of the peak amplitudes. It has long been recognized that the loca- 
tions and amplitudes of the peaks in the region 90 < £ < 900 (where the 
anisotropies are related to causal processes occurring in the photon-baryon 
plasma until recombination) are very sensitive to the variations in the pa- 
rameters of the model and hence serve as a sensitive probe to constrain the 
cosmological parameters and discriminate among various models (Hu et al. 
2001; Doran & Lilley 2001, 2002). In fact, for £ > 40, the ratios of the peak 
amplitudes are insensitive to the intrinsic amplitude of the CMB spectrum. 
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This renders the positions of the peaks, particularly the position of the first 
peak, as a powerful probe of the parameters of the model. The first-year 
observations of WMAP have measured the position of the first peak very 
accurately at i = 220.1 ± 0.8 (1 a) (Page et al. 2003), which we shall use in 
our fit. 

The angle 9 a, which is given by the ratio of sound horizon to the distance 
(angular diameter distance) of the last scattering surface, sets the acoustic 
scale £a through 

tx_ it S n(z dec ) 

0a JZ c c s (z)dz/H(zy [b) 

where the speed of sound c s in the plasma is given by c s = 1/ ^3(1 + R) and 
R = 3pb/4p 7 = 3f2bo/[4^ 7 o(l + z )] corresponds to the ratio of baryon density 
to photon density. The location of i-th peak in the angular power spectrum 
is given by 

^peaki = Mi - ( 16 ) 

where the phase shift Si, caused by the plasma driving effect, is determined 
predominantly by the pre- recombination physics (Hu et al. 2001) and can 
be approximated by 

*' = a '| o.3 n m0 } • (17) 

where a« is respectively 0.267, 0.24 and 0.35 for peak 1; peak 2 and peak 3 . 
Let us recall that f2 r0 gets contributions from photons (CMB) as well as 
from neutrinos, i.e., Q r = f2 7 + Q u . The present photon contribution to the 
radiation can be estimated from the CMB temperature T = 2.728K. This 
gives f2 7 o ~ 2.48 hr 2 x 10~ 5 , where h is the present value of the Hubble 
parameter in units of 100 km s _1 Mpc -1 . The neutrino contribution follows 
from the assumption of 3 neutrino species, a standard thermal history and a 
negligible mass compared to its temperature (Hu & Dodelson 2002) leading 
to tt^ w 1.7 h- 2 x 10- 5 . 

Equations (15)-(17), suplimented by (8) and (9), are now fully capable 
to compute the locations of the peaks for given values of the free parameters 
fl m o, fl<po, &w and h. We notice that a sufficiently big range of these param- 
eters produce the £ P eaki values in the observed range. For example, by fixing 
f2t,o — 0.05 and h = 0.65, the following choices of fl m0 and fl^o yield: 

Vl m0 = 0.3, = 0.6 -> £ peakl = 222.5, £ pC ak 2 = 536.7, £ pC ak 3 = 808.2; 
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n m0 = 0.25, fy = 0.65 -> ! pcakl = 225.0, £ pC ak 2 = 545.1, £ pC ak 3 = 820.9; 

Q m0 = 0.2, ^ = 0.72 -> £ peakl = 220.0, £ peaka = 540.7, £ pC ak 3 = 814.3, 

which are in reasonable agreement 1 with £ pcakl = 220.1 ± 0.8 and £ pcak2 — 
546 ± 10 measured by WMAP. These measurements are also consistent with 
other observations; for example, the BOOMERANG project (de Bernardis 
et al., 2002) measured the first three peaks in the ranges: 
£ peakl : 200 - 223, £ pcak2 : 509 - 561, £ pcak3 : 820 - 857 at 68% confidence 
level; and 

£ P ea kl : 183-223, £ pC ak 2 : 445-578, £ pC ak 3 : 750-879 at 95% confidence level, 
which are also in agreement with many other observations like, MAXIMA, 
DASI, CBI, etc., at least on the location of the first peak (Lee et al. 2001; 
Halverson et al. 2002; Sievers et al. 2002). 

In Figure 2, we have shown the allowed region in the Q m o — Q^o plane 
which produces the first peak £ peakl — 220.1±0.8 (68% confidence level). The 
marginalization over the other parameters Q^q and h is achived by taking 
projection of the full 4-dimensional region on the f2 m0 — f^o plane (Press et 
al 1986). 

It may be noted that the WMAP results obtained by Spergel et al.(2003)- 
that only < —0.78 can be accommodated within 95% confidence regions 
- comes from a lot of assumptions, some of which are not consistent with 
many observations. For example, there are several observations which also 
measure smaller values of H , apart from the higher values (see section 4.4). 
However, this degeneracy has not been taken into account and they consider 
only that HST observation which gives H = 0.72 ± 3(stat)±7(systematic) 
km s _1 Mpc -1 (Friedman et al. 2001) close to their best fit value. Note 
that there is also another HST Key Project which gives H = 64tg km s _1 
Mpc -1 (Saurabh et al. 1999). Sandage and his collaborators find a value 
even as low as Ho = 58 ± 6 km s" 1 Mpc -1 from an analysis of SNe la 
distances (Parodi et al. 2000). In Figure 2, we have instead considered all 
those possible combinations of the parameters f2 m0 e [0 — 1], f^o G [0 — 1.5], 
tt h0 E [0 - 0.1] and h E [0.5 - 0.8, which produce £ pcakl = 220.1 ± 0.8. This 

^^One can have even better fit by adjusting the parameters finely. For example, one can 
have 

r! m0 = 0.3, r> 0o = 0.61 £ peakl - 219.8, £ pcak2 = 530.2, ^ pcak3 = 798.4; 
O m0 = 0.25, S1 0O = 0.666 -» £ pcakl = 220.3, £ pcak2 = 533.6, £p eak3 = 803.5; 
with the same fibo and h. 
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amounts to a range of flwh 2 as [0 — 0.06] which safely contains the range 
0.01 < Qboh 2 < 0.04 which one essentially requires to explain the observed 
abundances of helium, deuterium and lithium (Narlikar & Padmanabhan 
2001). It should also be noted that the most favoured values like f2 m0 ~ 0.25 
and fibo^ 2 ~ 0.02 (or Q^ m 0.05) are compatible with a moderate h ph 0.65 
in this model, as we have shown in our examples. 




Figure 2: The blue region shows the allowed ranges of O m o and Q^q by the WMAP 
which produce the first peak in the range 220.1 ±0.8. This has been obtained after 
marginalizing over fibo and h (which have been varied in the ranges [0 — 0.1] and 
[0.5 — 0.8] respectively). The red-shaded contour corresponds to 95% confidence 
region (marginalized over Ai) from the SNe la data (with all 57 points, including 
the two newly discovered points SN 2002dc and SN 2002dd) and the red curve 
corresponds to the boundary of the 99% confidence region given by the same data 
(in calculating these regions, the extinction of SNe light by metallic whiskers, 
as discussed in section 5, has not been taken into account). The line in black 
corresponds to the flat model J7 m o + fi^o = !• 
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4.3 SNe la Observations 

For the SNe la data in the present model, y 2 decreases for lower values of 
Q m0 and gives the minimum x 2 f° r negative Q m0 . For example, for the flat 
model, the minimum value of x 2 is obtained as 62.93 for f2 m0 = —0.31 at 
53 dof (from the older sample of 55 SNe) and 63.47 for f2 m0 = —0.30 at 55 
dof (with the addition of the new points), which are though not physical. 
A 'physically viable' best-fitting solution, in this case, can be regarded as 
tt m0 = 0, which gives xV dof = 64.56/53 = 1.22 for = 0.14 and M = 
24.03 with Q = 13.3%, from the older sample, which represents a reasonably 
good fit, though not as good as the best fit in the standard cosmology with a 
constant A. Moreover, the addition of the new points to this sample improves 
the fit, giving x 2 /dof = 65.02/55 = 1.18 for ^ = 0.23 and M = 24.04 with 
Q = 16.7%, as the new best-fitting solution. 

The allowed region by the data is sufficiently large, as shown in Figure 
2, and low density models, which are also consistent with the WMAP ob- 
servation, are easily accommodated within the 95% confidence region (Note 
that the plotted allowed region by WMAP is only at la level. At higher 
levels, the region will be wider). For example, the following models represent 
reasonable fit. 

Ono = 0.3, {1^ = 0.6: x 2 /dof = 71.93/54 = 1.33 with Q = 5.2%; 

tt m0 = 0.25 = 1 - IV X 2 /dof = 71.36/54 = 1.32 with Q = 5.7%; 

n m0 = 0.25, fy = 0.65: x 2 /dof = 70.76/54 = 1.31 with Q =6.3%; 

tt m0 = 0.2, tt^o = 0.7: x 2 /dof = 69.63/54 = 1.29 with Q = 7.5%, etc., which 

are obtained from the older sample of 55 SNe. Addition of the new points 

improves the fits, giving: 

Ono = 0.3, IV, = 0.6: xVdof = 72.44/56 = 1.29 with Q = 6.9%; 
tt m0 = 0.25 = 1 - X 2 /dof = 71.87/56 = 1.28 with Q = 7.5%; 
^ m0 = 0.25, tl^o = 0.65: x 2 /dof = 71.22/56 = 1.27 with Q =8.3%; 
tt m0 = 0.2, tt^ = 0.7: x7 dof = 70.06/56 = 1.25 with Q = 9.8%, etc. 
One can go up to even as high as f2 m0 ~ 0.4 at 99% confidence level. 

4.4 Age of the Universe 

The parameters H and fl m0 set the age of the universe in this model. Re- 
member that there is still quite large uncertainty in the present value of Hq. 
Sandage and his collaborators find Hq = 58 ± 6 km s^ 1 Mpc -1 from an anal- 
ysis of SNe la distances (Parodi et al. 2000). This is also consistent with the 
value obtained from an analysis of clusters using Sunyaev-Zeldovich effect 
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which gives H = 60 ± 10 km s" 1 Mpc' 1 (Birkinshaw 1999). An HST Key 
Project supplies H = 64+6 km s" 1 Mpc -1 (Saurabh et al. 1999). Some 
experiments also measure higher H , for example, another HST Key Project, 
which uses Cepheids to calibrate several different secondary distance indica- 
tors, finds H = 72 ± 3(stat)±7(systematic) km s _1 Mpc -1 (Friedman et al. 
2001). Also by using the fluctuations in the surface brightness of galaxies 
as distance measure, it has been found that H = 74 ± 4 km s _1 Mpc -1 
(Blakeslee et al. 1999). Thus it appears that H lies somewhere in the range 
(50 — 78) km s _1 Mpc -1 . An average value of Hq = 65 km s _1 Mpc -1 from 
this range, constrains fl m o of the model by < Q m o < 0.34 to give the age 
of the universe to > 12 Gyr, so that the age of the oldest objects detected 
so far, e.g., the globular clusters of age tec = 12.5 ±1.2 Gyr (Cayrel et 
al. 2001; Gnedin et al. 2001), can be explained. The presently favoured 
value Q m0 ^0.25 easily accommodates in this range. It is interesting to note 
that the average value of H Q = 65 km s" 1 Mpc -1 we have preferred, is in 
good agreement with Hq = 67 ± 2(stat)±5(systematic) km s _1 Mpc^ 1 ob- 
tained recently by Gott et al (2001) from an analysis based on the median 
statistics. 

5 EXTINCTION BY METALLIC DUST 

In this section, we shall discuss the absorption of light by metallic dust 
ejected from the SNe explosions — an issue which is generally avoided while 
discussing m-z relation for SNe la. Although a number of observers believe 
this to be not a significant effect (see, for example, Riess et al (2001)), how- 
ever, taking this effect into consideration does improve the fit to the data, as 
we shall see in the following. 

It is well known that the metallic vapours are ejected from the SNe ex- 
plosions which are subsequently pushed out of the galaxy through pressure 
of shock waves (Hoyle & Wickramasinghe, 1988; Narlikar et al, 1997). Ex- 
periments show that metallic vapours on cooling, condense into elongated 
whiskers of ^0.5 — 1 mm length and ^10 -6 cm cross-sectional radius (Hoyle 
et al, 2000). Indeed this type of dust extinguishes radiation travelling over 
long distances (Aguire, 1999; Vishwakarma, 2002a). The density of the dust 
can be estimated along the lines of Hoyle et al (2000). If the metallic whisker 
production is taken as 0.1 M per SN and if the SN production rate is taken 
as 1 per 30 years per galaxy, the total production per galaxy (of spatial den- 
sity w 1 per 10 75 cm 3 ) in 10 10 years is ^ 2/3 x 10 41 g. The expected whisker 
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density, hence, becomes 2/3 x 10 41 x 10~ 75 ~ 10~ 34 g cm" 3 . We shall later 
see that this value is in striking agreement with the best-fitting value coming 
from the SNe la data. 

In an isotropic and homogeneous universe, the contribution to the ef- 
fective magnitude arising from the absorption of light by the intervening 
whisker-like dust, is given by 

rl(z) rz q}^' 

Am(z) = J o K Pg dl = K Pg0 J o (1 + z'Yj^y (18) 

where k is the mass absorption coefficient, which is effectively constant over 
a wide range of wavelengths and is of the order 10 5 cm 2 g _1 (Wickramasinghe 
& Wallis, 1996), p g oc S* -3 is the whisker grain density and l(z) is the proper 
distance traversed by light through the inter-galactic medium emitted at the 
epoch of redshift z. The net magnitude is then given by 

m QCt (z) = 771(2) + Am(z), (19) 

where the first term on the r.h.s. corresponds to the usual magnitude from 
the cosmological evolution given by equation (10). 

We note that taking account of this effect improves the fit to the SNe la 
data considerably. For example, the model with f2 m0 = 0.25, Q^o = 0.65 now 
gives x 2 /dof = 65.25/53 = 1.23 for M = 24.03 and p g0 = 2.02 x 10~ 34 g 
cm -3 with Q = 12.1%, from the older sample of 55 points. The addition of 
the new points improves the fit further by giving x 2 /dof = 65.75/55 = 1.20 
for M = 24.03 and p g0 = 1.89 x 10" 34 g cm" 3 with Q = 15.2%. Even 
the Einstein-deSitter model (f2 m0 = 1, fi^o = 0) gives an acceptable fit: 
x 2 /dof = 68.40/53 = 1.29 for M = 24.04 and p g0 = 4.95 x 10" 34 g cm" 3 
with Q = 7.6% (from the 55 points-data) and x 2 /dof = 68.97/55 = 1.25 for 
M = 24.04 and p g0 = 4.75 x 10" 34 g cm" 3 with Q = 9.8% (by adding the 
new points). Interestingly, this model (Einstein-deSitter) is also consistent 
with the CMB observations: Q h0 = 0.05 and h = 0.65 yield £ pcakl = 202.9, 
^peak 2 — 476.8, £peak 3 = 717.6. One can improve the fit by increasing fl h0 
and/or decreasing h: i\o = 0.1 and h = 0.55 yield £ pea ki = 220.4, £ P eak 2 — 
521.4, 

Speaks — 784.9. Also a better fit can be obtained in open models. For 
example, the model Q ni0 = 0.87, f2^ — with fl b0 = 0.05 and h = 0.65 
yields £ pea ki = 220.1, ^ P cak 2 — 518.8, £ pca k 3 = 780.9. This model also has an 
acceptable fit to the SNe la data: xV dof = 67.89/53 = 1.28 for M = 24.04 
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and p g o = 4.34 x 10 34 g cm 3 with Q = 8.2% (from the 55 points-data) and 
X 2 /dof = 68.45/55 = 1.24 for M = 24.04 and p g0 = 4.15 x 10~ 34 g cm" 3 with 
Q = 10.5% (by adding the new points). However, these models suffer from 
the age problem if h is not sufficiently low. For example, h should be < 54 
for Q m0 = 1 with Q^q = 0. Additionally, there is much evidence for low fi m o, 
as reviewed by Peebles & Ratra (2003). 

6 CONCLUSION 

There seems to be an impression in the community that the current obser- 
vations, particularly the high redshift SNe la observations and the measure- 
ments of the angular power fluctuations of the CMB, can be explained only 
in the framework of an accelerating universe. This, however, does not seem 
correct. The allowed parameter space by the datasets is wide enough to 
accommodate decelerating models also. We have shown that both these ob- 
servations can also be explained in a decelerating low density-model with a 
dark energy equation-of-state = —1/3 and the preferred curvature of the 
spatial section is slightly negative. For this equation of state, the resulting 
'dark energy' does not contribute to the expansion dynamics of the model 
(described by the Hubble parameter) and contributes to the curvature only. 

In order to fit the model to the data, we have considered the most recent 
observations. For example, for the SNe data, we have considered the older 
sample of 55 SNe of type la (54 SNe used by Perlmutter et al + SN 1997ff at 
z = 1.755) together with the two newly discovered supernovae SN 2002dc at 
z = 0.475 and SN 2002dd at z = 0.95. Addition of these new points to the 
older sample improves the fit to, more or less, all the models. For CMB, we 
have considered the first- year observations of WMAP which have measured 
the position of the first peak very accurately. 

It may be noted that the case = —1/3 is not special in any sense to the 
datasets and slightly more decelerating or slightly less decelerating models 
are also consistent with both the observations. In order to verify this, we 
change slightly from = —1/3 (in both directions) to, say, = —0.3 
and = —0.35 and check the status of the resulting models. The result is 
the following. 

= —0.3: 

A test model, for example, fl m0 = 0.25, f^o = 0.65, f2 b0 = 0.05 and h = 0.65 
yields £ pC aki = 220.2, £ P eak 2 — 533.4, £p ea k 3 = 803.2. The SNe la data, for the 
same model, give x 2 /dof = 72.10/54 = 1.34 with Q = 5.0%, from the older 
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sample of 55 points and x 2 /dof = 72.60/56 = 1.30 with Q = 6.7%, by adding 
the new points to this sample. 

= —0.35: 

The same test model, in this case, yields £ pca ki = 227.3, £ pca k 2 = 550.5, 
Speaks = 829.1. The SNe la data, for this case, give xV dof = 70.12/54 = 1.3 
with Q = 6.9% , from the older sample of 55 points and x 2 /dof = 72.57/56 = 
1.26 with Q = 9.1%, by adding the new points to this sample. These fits 
are acceptable and comparable to their respective values for = —1/3 
mentioned in sections 4.2 and 4.3. 

We also note that if we take into account the extinction of SNe light 
by the inter-galactic metallic dust, then the observed dimming of the high 
redshift SNe la can also be explained by the models without any dark energy, 
such as, the Einstein-deSitter model. These models are also consistent with 
CMB observations. In fact, there is a degeneracy in the fi m0 — f2^ plane 
along a line Q m0 + Q^q 1 and a wide range of Q m0 is consistent with the 
WMAP observation. 

Interestingly, another alternative explanation of the observed faintness 
of SNe la at large distances can be given in terms of a quantum mechan- 
ical oscillation between the photon field and a hypothetical axion field in 
the presence of extra-galactic magnetic fields. To satisfy other cosmological 
constraints, one then simply needs some form of uniform dark energy with 
Wfj) ~ —1/3 and the universe would be decelerating (Csaki et al. 2002). 

We conclude that it is premature to claim, on the basis of the existing 
data, that the present expansion of the universe is accelerating (or deceler- 
ating). Only more accurate SNe la data with z significantly > 1 can remove 
this ambiguity, as q is sensitive significantly to the SNe la data only. Whereas 
the CMB observations are consistent with both - accelerating as well as decel- 
erating models, as mentioned above. This endeavour may be accomplished 
by the proposed SuperNova Acceleration Probe (SNAP) experiment which 
aims to give accurate luminosity distances of type la SNe up to z ~ 1.7. 
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